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Superfluidity, the ability of a liquid or gas to flow with zero viscosity, is one of the most remarkable 
implications of collective quantum coherence. In equilibrium systems like liquid ^He and ultracold 
atomic gases, superfluid behaviour conjugates diverse yet related phenomena, such as persistency 
of metastable flow in multiply connected geometries and the existence of a critical velocity for 
frictionless flow when hitting a static defect. The link between these different aspects of superfluid 
behaviour is far less clear in driven-dissipative systems displaying collective coherence, such as 
microcavity polaritons, which raises important questions about their concurrency. With a joint 
theoretical and experimental study, we show that the scenario is particularly rich for polaritons 
driven in a three-fluid collective coherent regime so-called optical parametric oscillator. On the one 
hand, the spontaneous macroscopic coherence following the phase locking of the signal and idler 
fluids has been shown to be responsible for their simultaneous quantized flow metastability. On the 
other hand, we show here that pump, signal and idler have distinct responses when hitting a static 
defect; while the signal displays hardly appreciable modulations, the ones appearing in pump and 
idler are determined by their mutual coupling due to nonlinear and parametric processes. 

PACS numbers: 03.75.Kk, 71.36.+C, 42.65.Yj 


I. INTRODUCTION 

Microcavity polaritons, the novel quasiparticles result¬ 
ing from the coherent strong coupling between quantum 
well excitons and cavity photon^, have unique mixed 
matter-light properties that none of their constituents 
displays on its own. Because of their energy dispersion 
and their strong non-linearity inherited from the exci- 
tonic components, polaritons continuously injected by an 
external laser into a pump state with suitable wavevec- 
tor and energy can undergo coherent stimulated scatter¬ 
ing into two conjugate state^^J®, the signal and the idler, 
a process known as optical parametric oscillator (OPO). 
Since their first realisatiorP^, the interest in microcavity 
optical parametric phenomena has involved several fields 
of fundamental and applicative researclP^®^. 

Recently, considerable resources have been invested in 
exploring the fundamental properties of parametric pro¬ 
cesses, including the possibility of macroscopic phase co¬ 
herence and superfluid behavioui^. In spite of the coher¬ 
ent nature of the driving laser pump, the OPO process 
belongs to the class of non-equilibrium phase transitions, 
in which a 1/(1) phase symmetry is spontaneously bro¬ 
ker?^. While the phase of the pumped mode is locked 
to the incident laser, the ones of signal and idler are free 
to be simultaneously rotated in opposite directions. Be¬ 
cause of this phase freedom, recent experimenthave 
tested the OPO superfluid properties by exploring the 
physics of the signal-idler order parameter, demonstrat¬ 


ing the existence and metastability of vortex configura¬ 
tions. As the order parameter involves both signal and 
idler, their phase winding have opposite signs^^HSI] Qyu- 
cially, this makes both OPO fluids to display quantized 
flow metastability simultaneously. 

While in equilibrium condensates different aspects of 
superfluidity are typically closely related^, this is no 
longer true in a non-equilibrium context as for micro¬ 
cavity polariton^. In particular, those aspects of su¬ 
perfluidity related to the frictionless flow around defects 
are expected to be much more involved in OPO conden¬ 
sates than for any other investigated polariton co nden- 
sates, such as for the case of incoherent pumpin^^^^, 
and single-state resonantly pumped micro cavities^. In¬ 
dependently on the pumping scheme, the driving and 
the polariton finite lifetime prompt questions about the 
meaning of superfluid behaviour, when the spectrum of 
collective excitations is complex rather than real, rais¬ 
ing conceptual interrogatives about the applicability of a 
Landau criteriorP^. Yet, an additional complexity char¬ 
acterises the OPO regime, where the simultaneous pres¬ 
ence of three oscillation frequencies and momenta for 
pump, signal and idler correspondingly multiplies the 
number of collective excitation branche^^^. Note that 
from the experimental point of view, pioneering experi¬ 
ment have observed a ballistic non-spreading propaga¬ 
tion of signal/idler polariton wavepackets in a triggered- 
OPO configuration. However, given the complexity of the 
dynamics as well as the nonlinear interactions involved 
in this time-dependent conflguratiorP^, theoretical un- 



2 


derstanding of these observations is not complete yet. 

This Article reports a joint theoretical and experimen¬ 
tal study of an OPO configuration where a wide and 
steady-state condensate hits a stationary localised defect 
in the microcavity. Contrary to the criterion for quan¬ 
tized flow metastability for which signal and idler display 
simultaneous locked responses, we find that their scatter¬ 
ing properties when the OPO hits a static defect are dif¬ 
ferent. In particular we investigate the scattering prop¬ 
erties of all three fluids, pump, signal and idler, in both 
real and momentum space. We find that the modula¬ 
tions generated by the defect in each fluid are determined 
not only by its associated Rayleigh scattering ring, but 
each component displays additional rings because of the 
cross-talk with the other components imposed by non¬ 
linear and parametric processes. We single out three 
factors determining which one of these rings influences 
the most each fluid response: the coupling strength be¬ 
tween the three OPO states, the resonance of the ring 
with the blue-shifted lower polariton dispersion, the val¬ 
ues of each fluid group velocity and lifetime together es¬ 
tablishing how far each modulation can propagate from 
the defect. The concurrence of these effects implies that 
the idler strongly scatters inheriting the same modula¬ 
tions as the pump, while the modulations due to its own 
ring can propagate only very close to the defect and can¬ 
not be appreciated. Yet, the modulations in the signal 
are strongly suppressed, and not at all visible in exper¬ 
iments, because the slope of the polariton dispersion in 
its low momentum component brings all Rayleigh rings 
coming from pump and idler out of resonance. 

Note that the kinematic conditions for OPO are in¬ 
compatible with the pump and idler being in the sub¬ 
sonic regime. Thus, the coupling between the three com¬ 
ponents always implies some degree of scattering in the 
signal. In practice, the small value of the signal momen¬ 
tum strongly suppresses its actually visible modulations, 
something confirmed by the experimental observations. 


II. MODEL 

The dynamics of polaritons in the OPO regime, and 
their hydrodynamic properties when scattering against a 
defect, can be described via a classical driven-dissipative 
non-linear Gross-Pitaevskii equation (GPE) for the c on- 
pled exciton and cavity fields '0x,c {r,t) (^= 

The dispersive X- and C-fields decay at a rate 7x,c and 
are coupled by the Rabi splitting Qji, while the non¬ 
linearity is regulated by the exciton coupling strength 
9x- 
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FIG. 1. (Color online) OPO mean-field blue-shifts and fluctu¬ 
ation Rayleigh rings in the linear response scheme for homoge¬ 
neous pumping. Left panel: Signal s ([blue] upper triangles), 
pump p ([red] circles), and idler i ([green] lower triangles) 
mean-field energy blue-shifts en=s,p,i (in units of jp = 7 kp) 
versus the rescaled pump intensity Ip (in units of 7 ^) in the 
optical limiter regime. Parameters are Hr — 5 meV, zero 
cavity-exciton detuning, 7 x — — 0.12 meV, ujp — — 

— 1.25 meV, kp — 1.6/i?Ti“^, ks — 0, and ki — 3.2/im“^. The 
shaded area is stable OPO region, while the vertical dashed 
line corresponds to the pump power value chosen for plot ting 
the right panel. Right panel: Blue-shifted LP dispersion |To| 
with superimposed Rayleigh curves F^ ^ k+kp ^ evaluated 
within the linear response approximation (same symbols as 
left panel). The two rings corresponding to the signal state, 
F^ (u ^;) k 5 S'!*® shrinked to zero because ks — 0. 


We describe the defect via a potential V^(r) acting on the 
photonic component; this can eithe r be a d efect in the 
cavity mirror or a localised laser fiel d^^ * ^^ ^. In the con¬ 
servative, homogeneous, and linear regime (jx,c = 0 = 
Vrf(r) = gx)^ the eigenvalues of H are given by the lower 
(LP) and upper polariton (UP) energies, = 

T \J (<^k ~ cavity is driven by 

a continuous-wave laser field Fp(r, t) = 
into the OPO regime: Here, polaritons are continuously 
injected into the pump state with frequency ujp and 
momentum and, above a pump strength threshold, 
undergo coherent stimulated scattering into the signal 
{ujs^^s) and idler (co’i,k^) states. 

As a first step, it is useful to get insight into the system 
behaviour in the simple case of a homogeneous pump of 
strength J-p^) = fp. A numerical study of the coupled 
equations Q for the more realistic case of a finite-size 
top-hat pump profile .^^(r) will be presented later. To 
further simplify our analysis, we assume here that the UP 
dispersion does not get populated by parametric scatter¬ 
ing processes and thus, by means of the Hopfield coefh- 

dents 2X^, 2C^ = 1 ± - <^k )/(T “ W 

we project the GPE § onto the LP componentP^U 
X = XkV’x,k + C'kV’c.k, where t) = X^k 
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Here, 7k = is the effective LP decay 

rate, the interaction strength is given by 5 'k,ki,k 2 = 
^x^k^ki+k 2 -k^ki^k 2 , and the pumping term by 
fp{t) = Ck,/pe--^b 


III. LINEAR RESPONSE THEORY 

In the limit where the homogeneously pumped system 
is only weakly perturbed by the external potential Vrf(r), 
we apply a linear response analysi^^ The LP field is 
expanded around the mean-field terms for the three n = 
1, 2, 3 = s,p, i OPO stated, 

3 

n=l 

where k = k — k^. Eq. ^ is expanded linearly in both 
the fluctuations terms, u ~ and v as well as the de- 

n,k n,k 

feet potential. At zero-th order, the three complex uni¬ 
form mean-field equations can be solved to obtain the 
dependence of the signal, pump and idler energy blue- 
shifts, Cn = gxX^ the system parameter d^^ l ^^ l 

A typical behaviour of as a function of the rescaled 
pump intensity Ip = QxC^^fp in the optical limiter 
regime is plotted in the left panel of Fig.[^ At first order, 
one obtains six coupled equations diagonal in momentum 
spac^i^ 


^Jk,0 


■ U rC 
n,k 


■ V re 
n, —k 


iojt 
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= £kWk +\^d, (5) 

for the 6-component vector Wk = 
and for the potential part, = 

^,Ck^Ck+k.Ed(k), -Ca.a.-kHrf(-k))^. In 

^ we have only kept the terms oscillating at the 
frequencies ± cj and neglected the other terms in the 
expansion (i.e., 2ujn — oj^n i which are oscillating 
at frequencies far from the LP band, and thus with 
negligible amplitudes. In the particle-like and the 
hole-like channels, the Bogoliubov matrix determining 
the spectrum of excitations can be written a^^^ 


>Ck = 
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where the 3 OPO states components are 
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3 
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In absence of a defect potential = 0), Eq. ^ is the 
eigenvalue equation for the spectrum of excitations of a 
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FIG. 2. Linear responses to a static defect of the three OPO 
states in real and momentum space. Rescaled filtered OPO 
emissions (signal [top panels], pump [middle], and idler [bot¬ 
tom]) in real space |'0(r, ct;n)|^/|'0n|^ (left panels in linear 
scale) and momentum space (right panel in loga¬ 

rithmic scale) obtained within the linear response approxi¬ 
mation. The parameters are the same ones of Fig.[^ For the 
top left panel of the signal emission in real space, Gaussian 
filtering is applied to enhance the short wavelength modu¬ 
lations, which amplitudes are otherwise roughly 1% of the 
average signal intensity and about a factor of 10 times weaker 
than the modulation amplitudes in the pump fluid. 


homogeneous OPO, i.e., det(£^ — cj) = 0. The spectrum 
has 6 branches, uj^ (n v) k’ labeled by n = s,p, i and {u, v). 
Even though these degrees of freedom are mixed together, 
at large momenta, one recovers the LP dispersions shifted 
by the three states energies and momenta, i.e.,: 

- lim_ ^n,iu,v)M = M^k-kr. - ^n) , (9) 

k:^V‘^mc^R 

where + (—) corresponds to the u (v) particle- (hole- 
)hke branch. The OPO solution is stable (shaded area in 
Fig- 13 as far as k < 0- 

The shape of the patterns, or Cherenkov-like waves, 
resulting from the elastic scattering of the OPO 3-fluids 
against the static {uj = 0) defect can be determined start¬ 
ing from the spectrum, and in particular evaluating the 
closed curves F^ ^ in k-space, or “Rayleigh rings”!^ 

defined by condition (u v) k ~ "^he modu- 
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lations propagate with a direction fj^ ^ orthogonal 
to each curve a pattern wavelength given by 

the corresponding |k|, and a group velocity ~ = 

where = |v2„,„).k/^‘^n,(«,«),kl 

determines the distance, at any given direction f}^ k’ 
over which the perturbation extends away from the de¬ 
fect. For a single fluid under a coherent pump, the qual¬ 
itative shape of the modulation pattern generated in the 
fluid by t he de fect is mostly determined by the excitation 
spectrun j^^ l ^^ l 


For OPO, the spectrum of excitation on top of each of 
the three, n = 1, 2, 3, states (se^^ generates six identical 
Rayleigh rings F^ ^ for the 3 states. The Rayleigh 
rings for the OPO conditions specified in Fig.[2are clearly 
visible in the right panels of Fig. where we plot the 
k-space photoluminescence filtered at each state energy, 
i.e., \^|)iSMJn)? = l^/’n'^k.o + ^„,k + -kl^- 

chosen a ^-like defect potential, V^(k) = but we have 
however checked that our results do not depend on its 
exact shap^^. For the OPO conditions considered here, 
the signal momentum is at /Cg ^ 0, and thus only four 
of the six rings are present. The same rings are also 
plotted in the right panel of Fig. shifted at each of 
the three OPO state momentum k^, P^ k+k^ 
energies It is important to note that, even though 
the three OPO states have locked responses because the 
three states display the same spectrum of excitations, 
only one of the rings F^ is the most resonant 

at ujji with the interaction blue-shifted lower polariton 
dispersion. 


= uoi 


,LP 


+ 2X 


n=l 


( 10 ) 


where mean-field energy blue- 

shifts (measured in Fig. [T in units of 7^ = 7kp)- This 
implies that the most vi^le modulation for each fluid 
should be the most resonant one, with superimposed 
weaker modulations coming from the other two state 
rings. 

In the specific case of Fig.[^ the signal is at ks — 0 and 
thus produces no rings in momentum space. The other 
four rings are very far from being resonant with the blue- 
shifted LP dispersion (10) at cjs, and thus the signal dis¬ 


plays only an extremely weak modulation coming from 
the next closer ring, which is the one associated with the 
pump state, F^ ^ • We estimate that the signal mod¬ 

ulation amplitudes are roughly 1% of the average signal 
intensity and about a factor of 10 times weaker than the 
modulation amplitudes in the pump fluid. In order to 
show that the signal has weak modulations coming from 
the pump, we apply a Gaussian filter to the real space 
images (see inset of the left top panel of Fig. |^. As ex¬ 
plained in more details irP^, Gaussian filtering consists 
of subtracting from the original data a copy which has 
been convoluted with a Gaussian kernel, thus getting rid 


of the long-wavelength modulations. This procedure re¬ 
veals that indeed the pump imprints its modulation also 
into the signal, even though these are extremely weak, 
thus leaving the signal basically insensitive to the pres¬ 
ence of the defect. 

Pump and idler states are each mostly resonant with 
their own rings, i.e., q^^+k^ at ujp and T.^^+k, at 
cji, respectively. Thus one should then observe two su¬ 
perimposed modulations in both pump and idler filtered 
emissions, the stronger one for each being the most res¬ 
onant one. However, the modulations associated to the 
idler only propagate very close to the defect, at an av¬ 
erage distance ^ k before getting damped, 

and thus are not clearly visible. For the OPO conditions 
considered, this is due to the small idler group velocity 
~, as the dispersion is almost excitonic at the idler 

2,u,k 

energy. 

We can conclude that, for the typical OPO condition 
with a signal at ks — 0, considered in Figs. and the 
signal fluid does not show modulations and the extremely 
weak scattering inherited from the pump state can be 
appreciated only after a Gaussian filtering procedure of 
the image. In contrast, the idler has a locked response to 
the one of the pump state. Note that, for the conditions 
shown in Fig. as well as the other cases considered in 
Ref.l^, the subsonic to supersonic crossover of the pump- 
only stat^^ happens well above the region of stability of 
OPO. Thus it is not possible to study a case where the 
pump is already subsonic and at the same time promotes 
stimulated scattering. 


IV. EXPERIMENTS 

We now turn to the experimental analysis. We use a 
continuous-wave laser to drive a high quality {Q = 14000) 
GaAs microcavity sample into the OPO regime — de- 
tailson the sample can be found in a previous publica- 
tioi J^^PQ I. The polariton dispersion is characterised by a 
Rabi splitting Qr = 5.4 meV, the exciton energy ujq = 
1485.26 meV and we choose a sample region where the 
cavity-exciton detuning is slightly negative, —1 meV. We 
pump at kp = 0.89 /im“^ and ujp—uj^ = —2.43 meV, and, 
at pump powers 1.5-times above threshold, we obtain an 
OPO with signal at small wavevector kg = 0.21 /im“^ 
and ujg — uJq = —2.95 meV, and idler at ki = 1.57 /im“^ 
and uoi — ujq = —1.91 meV. The defect we use in the 
sample is a localized inhomogeneity naturally present in 
the cavity mirror. Note that the exact location of the 
defect can be extracted from the emission spectrum and 
is indicated with a dot (orange) symbol in the profiles of 

Fig.|3 

In order to filter the emission at the three states en¬ 
ergies, Is,p,i{T = x^y), and obtain 2D spatial maps for 
the OPO three states, we use a spectrometer and, at 
a fixed position xq, obtain the intensity emission as a 
function of energy and position, /(e,Xo,^). By chang¬ 
ing xq we build the full emission spectrum as a function 
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FIG. 3. (Color online) Experimental OPO spectrum and fil¬ 
tered emissions of signal, pump and idler in presence of a 
structural defect. The six panels show the filtered emission 
profiles in real space of signal (top), pump (middle) and idler 
(bottom) Is,p,i{r). A Gaussian hltering to enhance the short 
wavelength modulations is applied in the right column panels. 
Here, the extracted wave-crests from the idler emission (yel¬ 
low contours in the bottom panel) are also superimposed to 
the pump prohle (middle) by applying a 7r-phase shift. The 
(orange) dot indicates the position of the defect. The lower 
panel shows the experimental OPO spectrum. Energy and 
momentum of the three OPO states are labeled with a [blue] 
upper triangle (signal), a [red] circle (pump), and a [green] 
lower triangle (idler), while the localised state, clearly visi¬ 
ble just below the bottom of the LP dispersion, is indicated 
with the symbol d. The bare LP dispersion is extracted from 
an off-resonant low pump power measurement, as well as the 
emission of the exciton reservoir (X) and the one of the UP 
dispersion (each in a different scale). 


of energy and 2D position, /(e, r). The filtered emis¬ 
sion for each OPO state is obtained from the integrals 
In=s,p,i{^) = del{e^v)^ with cr = 0.08 meV. The 

results are shown in Fig. for respectively the signal 
(top panel), pump (middle) and the idler (bottom) pro¬ 


files. The signal profile shows no appreciable modulations 
around the defect locations, nor these could be observed 
after applying a Gaussian filtering procedure of the im¬ 
age. In contrast, in agreement with the theoretical re¬ 
sults, both filtered profiles of pump and idler show the 
same Cherenkov-like pattern. We extract the wave-crests 
from the idler profile ([yellow] contours in the bottom 
panel) and superimpose them to the pump profile (mid¬ 
dle panel) with an added 7r-phase-shift, revealing that the 
only modulations visible in the idler state are the ones 
coming from the pump state. 


V. NUMERICAL ANALYSIS 


The agreement between the results obtained experi¬ 
mentally and within the linear response approximation is 
additionally confirmed by an exact full numerical analy¬ 
sis of the coupled equations 0 for a finite size pump via 
a 5^^-order adaptive-step Runge-Kutta algorithm. De¬ 
tails are given iiP^. The pumping conditions are very 
similar to those previously considered in the linear re¬ 
sponse approximation of Figs, [^andj^ while the pump 
profile J>(r) is now finite-size top-hat. In particular, 
we consider the case of zero cavity-exciton detuning, 
kp = 1.6 /im“^. Up — u^ = —0.44 meV and the pump 
power strength is fixed just above threshold, so that to 
produce a stable steady state OPO with, in absence of 
the defect, signal at small wavevector kg = —0.2 /im“^ 
and idler at ki = 3.4 /im“^. 

When adding a localised defect potential, the steady 
state OPO develops Rayleigh rings in momentum space, 
yet, as shown irP^, the spectrum continues to be 6-\ike 
in energy, allowing to easily filter in energy the emis¬ 
sion of the three OPO states. Results are shown in 
Fig. where real-space emissions |'0c'(r, are plot¬ 

ted in the left panels, while the ones in momentum space 
|'0c'(k, ci;^)P in the right panels. We observe a very sim¬ 
ilar phenomenology to that one obtained in the linear 
approximation shown in Fig. The signal now is at 
slightly negative values of momenta kg = —0.2 /im“^, 
thus implying a very small Rayleigh ring associated with 
this state. Thus we observe that only the modulations 
associated with the pump are the ones that are weakly 
imprinted in the signal state and that can be observed by 
means of a Gaussian filtering (inset of top-left panel). We 
have fitted the upstream wave-crests and obtained the 
same modulation wave-vector as the pump one ([blue] 
upper triangles). Similarly to the linear response case, 
we also find here that the most visible perturbation in 
the emission filtered at the idler energy is the one due 
to the pump Rayleigh ring. As before, the modulations 
due to the idler Rayleigh ring cannot propagate far from 
the defect because of the small group velocity associated 
with this state. 
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FIG. 4. (Color online) Full numerical responses to a static de¬ 
fect of the three OPO states in real and momentum space. Fil¬ 
tered OPO emissions (signal [top panels], pump [middle], and 
idler [bottom]) in real space |'0c(r, Ct;n)|^ (left panels in linear 
scale) and momentum space |'0c(k, c^n)!^ (right panel in log¬ 
arithmic scale) obtained by a full numerical evaluation of 0 - 
For the top left panel of the signal space emission, Gaussian 
filtering is applied to enhance the short wavelength modula¬ 
tions of this state, revealing that the modulations correspond¬ 
ing to the pump state are also imprinted (though weakly) into 
the signal. The symbols indicate the pump ring diameter ex¬ 
tracted from fitting the upstream modulations and resulting 
in a density wave wave-vector coinciding with the one of the 
pump kp = 1.6 /im“^. 


the signal, because of its low characteristic wavevector, 
so much that experimentally cannot be resolved. The 
main features of the real- and momentum-space emission 
patterns are understood in terms of Rayleigh scattering 
rings for each component and a characteristic propaga¬ 
tion length from the defect; the rings are then trans¬ 
ferred to the other components by nonlinear and para¬ 
metric processes. Our theoretical and experimental re¬ 
sults stress the complexities and richness involved when 
looking for superfluid behaviours in non-equilibrium mul¬ 
ticomponent condensates such as the ones obtained in the 
optical parametric oscillation regime. 
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VI. CONCLUSIONS 


To conclude, we have reported a joint theoretical and 
experimental study of the superfluid properties of a non¬ 
equilibrium condensate of polaritons in the so-called op¬ 
tical parametric oscillator configuration by studying the 
scattering against a static defect. We have found that 
while the signal is basically free from modulations, pump 
and idler lock to the same response. We have highlighted 
the role of the coupling between the OPO components by 
non-linear and parametric processes. These are respon¬ 
sible for the transfer of the spatial modulations from one 
component to the other. This process is most visible in 
the clear spatial modulation pattern that is induced by 
the non-superfluid pump onto the idler, while the same 
modulations are only extremely weakly transferred into 
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In the manuscript, we carry on a theoretical and ex¬ 
perimental analysis of the response of microcavity polari- 
tons in the optical parametric oscillator (OPO) regime to 
a static defect. For the theoretical calculations we follow 
two independent approaches: In the first approach, we 
exactly numerically solve the dynamics of the two cou¬ 
pled Gross-Pitaevskii equations (GPEs) for the exciton 
and cavity fields for the realistic case of a finite-size top- 
hat pump profile In the second approach, we ap¬ 

ply a perturbative linear response approximation for the 
lower polariton (LP) state which leads to semi-analytical 
results in the limiting case of a spatially homogeneous 
pump profile. Both methods have been already success¬ 
fully used in the literature to explore several properties 
of the OPO operation. 

While all fundamental formulae and information has 
been given in the main text, in this Supplementary Ma¬ 
terial (SM) we present some additional details on both 
approaches that might be of interest to the specialized 
reader. In addition to that, we make use of the linear re¬ 
sponse approach to examine some regimes that are hardly 
accessed experimentally or within a full numerical ap¬ 
proach, which however allow to put the conclusions of 
our work into a broader perspective. 


VII. FULL NUMERICS 

The classical driven-dissipative non-linear Gross- 
Pitaevskii equation (GPE) for the coupled exciton and 
cavity fields ^x,c {r,t) {h=l) 


idt 


'ipx 


= H 


'Ipx 
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0 

Fp{r,t) 


( 1 ) 


where 


H 


+ gxli’xl 

Qr/2 


Q.r/2 




( 2 ) 

is solved numerically on a 2D grid oi N x N = 2^ x 2^ 
points and a separation of 0.47 /im (i.e., in a box L x 
L = 121 /imxl21 fim) by using a 5^^-order adaptive-step 
Runge-Kutta algorithm. Gonvergence has been checked 
both with respect the resolu tion i n space L/V as well as 
in momentum tt/L, withoufP^^ as well as in presence 
of the defect. The same approach has been already used 
in previous pub licatio ns from some of the authors (for a 
review, see Refs.^^^^^. As for the system parameters, we 
have considered a LP dispersion at zero photon-exciton 
detuning, uJq = uJq , a dispersionless excitonic spectrum, 
uj^ = uJq and a quadratic dispersion for photons = 
uJq /2mc^ with the photon mass me = 2.3 x 10“^me, 
where me is the bare electron mass. The LP dispersion. 


= uj^ + uj^ — y — 00 ^is characterised 
by a Rabi splitting Qji = 4.4 meV. Eurther, the exciton 
and cavity decay rates are fixed to yx = yc = 0.53 meV. 
Eor the defect we choose a (5-hke defect potential 

Vd(r) = gd6{r - vq) , (3) 

where its location fq is fixed at one of the N x N points 
of the grid. Note that in a finite-size OPO, local currents 
lead to inhomogeneous OPO profiles in within the pump 
spot, despite the external pump havi ng a top-hat pro¬ 
file with a completely flat inner regioiP®^ — as shown 
later, this can be observed in the filtered OPO profiles 
evaluated in absence of a defect shown as dashed lines 
in S. Eig. We have thus chosen the defect location so 
that it lies in the smoothest and most homogeneous part 
of the OPO profiles. 

Also we have checked that our results do not quali¬ 
tatively depend on the strength (nor on the sign) of 
the defect potential, as far as this does not exceed a crit¬ 
ical value above which it destabilises the OPO steady- 
state regime. The pump, Fp{r,t) = 
has a smoothen and rotationally symmetric top-hat pro¬ 
file, Fp{t) = Fp{r) = ^[tanh(f^^) -tanh(f^^)] with 
strength fp = 1.23/^^ = 0.053 meV/gm and param¬ 
eters ro = 8.68 /im, ap = 34.72 gm. We pump at 
kp = 1.6 gm~^ in the x-direction, = (A:p,0), and at 
Up — = —0.44 meV, i.e., roughly 0.5 meV above the 

bare LP dispersion. By increasing the pump strength 
/p, we find the threshold above which OPO switches 
on, leading to two conjugate signal and idler states. We 
then fix the pump strength just above this threshold 
{fp = 1.23/^^), where we find a steady state OPO so¬ 
lution which is stable (see Ref.l^ for further details). In 
absence of the defect, this condition corresponds to a sig¬ 
nal state at ks = —0.2 gm~^ and Us — Uq = —1.64 meV 
and an idler at ki = 3.4 gm~^ and Ui — = 0.76 meV. 

It is interesting to note that already very close to the 
lower pump power threshold for OPO, the selected sig¬ 
nal momentum is very close to zero. This contrasts with 
what one obtains in the linear approximation scheme, 
where instead just above the lower OPO threshold there 
exists a broad interval of permitted values for kg (and 
thus ki) — we will discuss further this “selection prob¬ 
lem” for parametric scattering later in the section “Linear 
response”. 

We evaluate the time dependent full numerical solution 
of 0 '0x,c(i*,t), until a steady state regime is reached. 
Here, both its Eourier transform to momentum k and 
energies uj can be evaluated numerically. We plot on 
the left panel of S. Eig. a cut at = 0 of the pho¬ 
tonic component of the OPO spectrum in presence of 
a point-like defect, |'0c'(/ca,, 0, ci;)p, as a function of the 
rescaled energy uj — uj^ versus the x-component of mo- 
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S. Fig. 1. (Color online) OPO spectrum obtained by full 
numerics. Left panel: Photonic component of the OPO spec¬ 
trum in presence of a point-like defect, |'0c(fc, 0, (loga¬ 
rithmic scale), as a function of the rescaled energy uj — ujq 
versus the x-component of momentum kx (cut at ky = 0) for 
a top-hat pump (see text for the space profile and parameter 
values), with intensity fp = 1.23/^^ above the OPO thresh¬ 
old, pump wave-vector kp = 1.6 in the a:-direction and 

^p — ^x — —0.44 meV. The symbols indicate the signal ([blue] 
upper triangle), pump ([red] circle), and idler ([green] lower 
triangle) energies, as well as the two momenta kx on each 
state Rayleigh ring at = 0. Note that the logarithmic 
scale results in a fictitious broadening in energy of the spec¬ 
trum, which is in reality (5-like (see right panel). The bare 
LP dispersion, including its broadening due to finite lifetime, 
is plotted as a shaded grey region, while the bare UP disper¬ 
sion as a (black) dot-dashed line. Right panel: Momentum 
integrated spectrum, |'0c(k, Ct;)|^ (linear scale) as a func¬ 
tion of the rescaled energy u — Uq , where it can be clearly 
appreciated that the emission is (5-like in energy. 


mentum kx (cut at ky = 0). In the right panel we plot 
instead the corresponding momentum integrated spec¬ 
trum, |'0c'(k, Cc;)p. Here, we can clearly see that the 
presence of the defect does not modify the fact that the 
OPO emission for the OPO signal ([blue] upper trian¬ 
gle), pump ([red] circle), and idler ([green] lower trian¬ 
gle) states has a completely flat dispersion in energy, thus 
indicating that a stable steady state OPO solution has 
been reached. Note that in the spectrum map of the left 
panel of Fig.[^ the logarithmic scale results in a fictitious 
broadening in energy. However, from the integrated spec¬ 
trum plotted in linear scale in the right panel of Fig.[2one 
can clearly appreciate that this emission is J-like, exactly 
as it happens for the homogeneous OPO casd^. Thus the 
effect of the defect is to induce only elastic (i.e., at the 
same energy) scattering; now the three OPO states emit 
each on its own Rayleigh ring (given each by the sym¬ 
bols on the left panel of S. Fig. which represent the 
rings at a cut for ky = 0). This makes it rather difficult 
to extract the separated signal, pump and idler profiles 
by filtering in momentum, as done previously for the ho¬ 
mogeneous case, but it still allows to filter those profiles 
very efficiently in energy. In fact, because the emission is 
(5-like, it is enough to fix a single value of the energy uj to 
the one of the three states 0Jn=s,p,i^ thus extracting the 
filtered profiles either in real space |'0c'(r, Cc;^)p or in mo¬ 
mentum space |'0c'(k, cjn)!^ — we have however checked 
that integrating in a narrow energy window around ujn 
does not quantitatively change the results. 
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S. Fig. 2. Real space signal, pump and idler OPO one¬ 
dimensional filtered profiles derived from finite size numer¬ 
ics with and without a defect. Rescaled OPO filtered emis¬ 
sions along the y — 6 direction, \'il;c{x^y — 0,ct;n)|^^^, ob¬ 
tained by numerically solving the GPE Eq. (1) of the main 
manuscript. While the dashed lines represent the filtered 
emissions of signal (top panel), pump (middle) and idler (bot¬ 
tom) for a top-hat pump without a defect, the solid lines are 
the same OPO conditions but now for a defect positioned at 
(xd, yd) = (9-5, —0.5) /xm corresponding to the vertical dotted 
lines. The system parameters are the same ones as those of 
Eig. 4 of the main manuscript. 


The results of the above described filtering are shown 
in Fig. 4 of the manuscript. In S. Fig. we instead 
plot the corresponding one-dimensional profiles in the 
^ = 0 direction both in presence (solid line) and with¬ 
out (dashed line) a defect. Here, we can observe that, 
even if the pump has a top-hat flat profile, as also com¬ 
mented previously, in absence of a defect, the finite-size 
OPO is characterised by inhomogeneous profiles of sig¬ 
nal, pump and idler because of localised currents. Fur¬ 
ther, we observe that the presence of a defect induces 
strong modulations in pump and idler. In order to re¬ 
veal that the pump also imprints its modulation into 





















the signal, even though these are extremely weak (and 
hardly visible in both the top main panel of Fig. 4 of 
the manuscript and the top panel of S. Fig. [^, we ap¬ 
ply a Gaussian filtering to the signal images, which re¬ 
sult is shown in the inset of the top-left panel of Fig. 4 
in the manuscript. This consists of the following proce¬ 
dure. The original data for the real space profile '0(r) are 
convoluted with a Gaussian kernel K{v — r'), obtaining 
a new profile, '0(r') = j dY%jj{Y)K{v — r'), where short 
wavelengths features are smoothen out. The convoluted 
image '0(r') is then subtracted from the original data, 
giving '0(r) — '^(r), and effectively filtering out all long 
wavelength details. The same procedure of Gaussian fil¬ 
tering is also applied to the signal emission profile ob¬ 
tained within the linear response analysis and shown in 
the inset of the top left panel of Fig. 2 of the manuscript. 



S. Fig. 3. Spectrum of collective excitations. Cut at 
/cy = 0 of the real part of the quasiparticle energy disper¬ 
sion plotted versus kx — kn- The spectrum 

is evaluated within the linear approximation scheme and the 
parameters are the same ones used for Fig. 2 of the main 
manuscript. 


VIII. LINEAR RESPONSE 

In the manuscript we also make use of the linear re¬ 
sponse approximation to analyse the OPO response to a 
static defect, valid for a homogeneous pump scheme. In 
order to do so, we first evaluate the Bogoliubov matrix 
given in Eq. (6) of the manuscript, whose eigenvalues de¬ 
termine the spectrum of collective excitations. We plot in 
S. Fig. 1^ a typical collective dispersion (here we consider 
the same system parameters as the ones used in Fig. 2 of 
the manuscript), by plotting the real part of the Bogoli¬ 
ubov matrix eigenvalues ^O0n,iu,v),\i-\in ^ function of 

kx — kfi (cut at ky = 0). Note that the Rayleigh rings can 
be found by finding the intersections 5Rco’^,(n,v),k-k^ = 0- 

As also done for the full numerics, we consider there a 
(5-like disorder potential Vrf(r) = gdS{Y — tq), and Figs. 1 
and 2 of the manuscript, as well as S. Figs. and be¬ 
low, are plotted for this case. We have however checked 
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S. Fig. 4. (Color online) OPO for finite and positive signal 
momentum in the linear response approximation for homo¬ 
geneous pumping. Rescaled OPO emission within the linear 
response approximation in real |'0(r,ct;n)|^/|'0n|^ (left panels 
in linear scale) and momentum space \^^{ujn)\^ (right panels 
in linear scale) filtered at the energies of the three OPO states: 
signal (top panels), pump (middle), and idler (bottom). The 
parameters are the same as those used in Fig. 2 of the main 
manuscript, with the exception of the signal and idler mo¬ 
menta, which are now at ks = 0.7 and ki = 2.5 

respectively. Each scale of variation for the profiles is plotted 
in the color-box next to the left panels. A cut in the y = 0 
direction for the three profiles is plotted in the bottom panel 
of S. Fig.[^ 


that our results do not depend on the specific shape of 
the defect potential: In particular we have also consid¬ 
ered the response to defects with smoothen Gaussian-like 
profiles, whose effect is only to partially weaken the up¬ 
stream modulations emission in real space. 

We have seen in the manuscript that both experimen¬ 
tally as well as in the full numerical analysis, one obtains 
OPO conditions where the signal momentum is very close 
to zero, kg ~ 0. The reason why OPO selects almost zero 
momentum signal conditions is still awaiting a theoreti¬ 
cal explaination. In contrast, the linear response approx¬ 
imation allows to access very different OPO mean-field 

































S. Fig. 5. (Color online) OPO for finite and negative signal 
momentum in the linear response approximation for homo¬ 
geneous pumping. Rescaled OPO emission within the linear 
response approximation in real |'0(r,ct;n)|^/|'0n|^ (left panels 
in linear scale) and momentum space |(c^n) | ^ (right panels 
in linear scale) hltered at the energies of the three OPO states: 
signal (top panels), pump (middle), and idler (bottom). The 
parameters are the same as those used in Fig. 2 of the main 
manuscript, with the exception of the signal and idler mo¬ 
menta, which are now at ks = —0.4 and ki = 3.6 /im“^ 

respectively. A cut in the y = 0 direction for the three prohles 
is plotted in the top panel of S. Fig. 


conditions, for which, by leaving unaltered the value of 
the pump momentum, the signal can appear at finite 
momentum values, either positive or negative. This is a 
peculiarity of this approximation scheme, where one can 
show that, at mean-field level, one has the possibility of 
choosing different values of the signal momentum (and 
thus the one of the idler k^), and that this choice range is 
quite broad particularly when the pump power is close to 
the lower pump threshold for OPO — the stability region 
for OPO is plotted as a shaded grey region in the top left 
panels of S. Figs. |^and[^ This is a well known “selection 
problem” for parametric scattering (see Ref.^^i The rea¬ 
son why in full numerical calculations, parametric scat¬ 
tering processes select a signal with a momentum close 
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S. Fig. 6. (Color online) Real space signal, pump and idler 
OPO one-dimensional hltered prohles in the linear response 
approximation. OPO hltered emissions along the y = 0 di¬ 
rection, I{x,y = 0,ct;n) = \'ip{x,y = 0,ct;n)|^/|'0n|^ rescaled by 
the mean-held solution in absence of the defect, for the 
three cases analysed above within the linear response approx¬ 
imation: the case of a signal at kg = —0.4 (top panel) 

corresponds to the same conditions as S. Fig. [^ , a signal at 
kg = 0.0 yin~^ (middle panel) corresponds to Fig. 2 of the 
main manuscript, and a signal at kg = 0.7 yin~^ (bottom 
panel) corresponds to S. Fig.[^ In each panel we plot, on the 
same linear scale, the hltered prohles of signal ([blue] dashed 
line), pump ([red] solid), and idler ([green] dotted), while the 
horizontal gray dashed lines represent the values of the mean- 
held emission prior to adding a defect (rescaled here to 1). 


to zero, already very close to the lower pump thresold for 
OPO, is still awaiting for an explaination. In particular, 
this quest cannot be addressed within a spatially homo¬ 
geneous approximation where the three mean-held solu¬ 
tions for pump, signal, and idler states can be describes 
by plane waves. One can only show that, within the 
same mean-held approximation scheme, when increasing 
the pump power towards the upper threshold for OPO, 
the blue-shift of the LP polariton dispersion due to the 








































S. Fig. 7. Real space profiles of three uncoupled signal, pump 
and idler fluids. The rescaled profiles |'?/;(r, ct;n)|^/|'0n|^ are ob¬ 
tained by setting all off-diagonal couplings in Eq. (6) of the 
main manuscript to zero, resulting in three uncoupled signal 
(left column), pump (middle column), and idler (right col¬ 
umn) fluids. The three rows correspond to the three different 
cases previously analysed within the linear response approx¬ 
imation: the case of a signal at kg = —0.4 /xm”^ (top row) 
corresponds to the same conditions as S. Fig. |^, a signal at 
ks = 0.0 /im“^ (middle row) corresponds to Fig. 2 of the 
main manuscript, and a signal at kg = 0.7 /xm“^ (bottom 
row) corresponds to S. Fig. 


increasing mean-field polariton densi^ causes the signal 
momentum to converge towards zercP^ ^ 0. 

OPO conditions with a finite momentum signal — 
Given the freedom of choice for the signal momentum 
close to the lower OPO threshold, we consider here two 
additional cases, that could not be studied neither exper¬ 
imentally, nor within a full numerical approach, but that 
instead we can easily analyse within the linear response 
theory. In particular, we have left fixed the pumping con¬ 
ditions kp = 1.6 /im“^ and uOp — ujx = —1-25 meV and 
considered two opposite situations. 

In the first case, the signal has a finite and positive 
momentum kg = 0.7 /im“^ and thus the idler is at low 
momentum, ki = 2.5 /xm”^. The results are shown in 
s. Fig. |§ Here we see that all six Rayleigh rings are 
clearly visible and, in addition, as the idler is at lower 
momentum compared to the case considered in Figs. 1 
and 2 of the manuscript, and thus its dispersion steeper, 
the idler group velocity is large enough to appreciate the 
modulation of the Rayleigh ring associated to this state. 
As a result, each of the three filtered OPO emissions 
exhibits as the strongest modulation the one coming from 
its own Rayleigh ring, included the signal which is now 
at finite and large momentum. In this case, the OPO 
response of each filtered state profile looks completely 
independent from the other, as if we were pumping each 


state independently. 

In the second case shown in S. Fig.|^, the signal is finite 
and negative, kg = —0.4 /xm“^, and the idler is now at 
very large momentum, ki = 3.6 /xm“^, where its disper¬ 
sion is very exciton-like and flat, and thus the idler has 
a very small group velocity and its own modulations are 
visible only very close to the defect. For this case, we can 
appreciate in the idler filtered profile overlapped modu¬ 
lations from all the three state Rayleigh rings (note that 
because the signal is at negative momentum, its modu¬ 
lations have an opposite direction compared to the ones 
of pump and idler), while in the signal we can mostly see 
the signal long wavelength modulations and only very 
weakly the pump one. 

We can compare the different modulation strengths 
of the three OPO profiles by looking at the color bars 
plotted next to the profiles. In order to better compare 
them on the same plot, we show in S. Fig. the one¬ 
dimensional OPO filtered emissions along the ^ = 0 di¬ 
rection for the three cases analysed above. While for the 
OPO conditions with a signal at kg = 0.0 /xm“^ (middle 
panel, corresponding to Fig. 2 of the main manuscript), 
the imprinted modulations from the pump are hardly vis¬ 
ible, and can only be appreciated after a Gaussian filter 
manipulation, both OPO cases with a finite momentum 
signal result in modulations in the signal with an ampli¬ 
tude of the same order of magnitude of both pump and 
idler fluids. 

In S. Fig. we instead show the real space profiles 
obtained by setting all off-diagonal couplings in the Bo- 
goliubov matrix of Eq. (6) of the main manuscript to 
zero, resulting in three uncoupled signal (left column), 
pump (middle column), and idler (right column) fluids. 
This underlines the importance of the coupling between 
the three fluids in the three different OPO regimes previ¬ 
ously analysed within the linear response approximation. 
In particular, for the OPO conditions such that the sig¬ 
nal has a finite and positive momentum k^ (bottom row 
of S. Fig. 7, which corresponds to the conditions shown 
in S. Fig. 4), the coupling has little effect and the three 
fluid respond to the defect in practice in a independent 
way (each Rayleigh ring influences its own fluid). The 
OPO condition with a finite and negative momentum k^, 
which corresponds to the conditions shown in S. Fig. 
are shown in the top row of S. Fig. Here, we can see 
that the coupling between the three fluids plays an effect. 
The modulations of the pump can be appreciated both in 
the signal profile (though weakly), as well as in the idler 
profile. In the idler fluid its own modulations can only 
propagate very close to the defect and thus the only re¬ 
ally appreciable modulations are the ones inherited from 
the pump. Finally, for the experimentally relevant case 
of an OPO with a signal at zero momentum (middle row 
of S. Fig. which corresponds to the conditions shown 
in Fig. 2 of the manucript), there is also a role played by 
the coupling, though, as already thoroghly analysed, the 
modulations in the signal inhereted from the pump can 
only be appreciated after a Gaussian filtering manipula- 













































tion. 

Finally note that, as it also happens for the OPO con¬ 
ditions shown in Fig. 1 of the manuscript, the subsonic 
to supersonic crossover of the pump-only stat^^ happens 
at pump intensities well above the region of stability of 


OPO — the shaded gray regions of S. Figs. andThus, 
it is not possible to study a case where the pump is al¬ 
ready subsonic and at the same time promotes stimulated 
scattering. 
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